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We consider a generalization of the Lemaitre-Tolman-Bondi (LTB) solutions by keeping the LTB 
metric but replacing its dust matter source by an imperfect fluid with anisotropic pressure Yl a b- 
Assuming that total matter-energy density p is the sum of a rest mass term, p*-" 1 ' , plus a radiation 
p^' = 3p density where p is the isotropic pressure, Einstein's equations are fully integrated without 
having to place any previous assumption on the form of Tl a b- Three particular cases of interest are 
contained: the usual LTB dust solutions (the dust limit), a class of FLRW cosmologies (the homoge- 
neous limit) and of the Vaydia solution (the vacuum limit). Initial conditions are provided in terms 
(iPsWable averages and contrast functions of the initial densities of p^ m ' , p' 1 "' and the 3-dimensional 
Ricci scalar along an arbitrary initial surface t = U. We consider the source of the models as an in- 
tSSctive radiation-matter mixture in local thermal equilibrium that must be consistent with causal 
Extended Irreversible Thermodynamics (hence Il a t is shear viscosity). Assuming near equilibrium 
conditions associated with small initial density and curvature contrasts, the evolution of the models 
isQbialitatively similar to that of adiabatic perturbations on a matter plus radiation FLRW back- 
gtmmd. We show that initial conditions exist that lead to thermodynamically consistent models, 
but only for the full transport equation of Extended Irreversible Thermodynamics. These interac- 
trW mixtures provide a reasonable approximation to a dissipative 'tight coupling' characteristic of 
radiation-matter mixtures in the radiative pre-decoupling era. 
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I. INTRODUCTION 



The Lemaitre-Tolman-Bondi (LTB) solutions with a dust source Jfij , M , 
H are widely popular models of cosmological inhomogeneities (see Igffbr a 
comprehensive review). We present a generalization in which the LTB met- 
ric is kept but the source is replaced by an imperfect fluid with anisotropic 
pressure, under the assumption that matter-energy density is decomposed 
Q-i as p = p( m ' + p^ r ' , a mixture of "rest mass" and "radiation" components 
(a mixture of non-relativistic and relativistic matter), so that = 3p, 
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where p is the isotropic pressure. The purpose of this paper is to derive 
the important geometric properties of the solutions within a convenient 
framework and then to examine the compatibility of the models with the 
k> , physics of radiation-matter sources. 

The study of inhomogeneous cosmological models such the LTB models 
is a well motivated and justified endeavor. First of all, it complements 
the usual perturbative approach by allowing one to study the non-linear 
evolution of inhomogeneities. Also, a nearly isotropic Cosmic Microwave 
Background Radiation (CMB) does not rule out an inhomogeneous universe 
compatible with current CMB observations [§]],§■ Moreover, the models 
we examine here are of particular interest when it is necessary to consider 
the radiation component and the dissipative processes associated with its 
interaction with non-relativistic matter 0, On the other hand, 

even if it is reasonable to use a dust source model as a theoretical matter 
model for present day universe, this source cannot describe the CMB and 
does not allow one to deal with temperatures of the photon gas. 

Various physical interpretations and observational bounds have been pro- 
posed for anisotropic stresses in cosmic matter sources [[To] , [jll| . If the mix- 
ture components are not interacting, then we have a mixture of collision-less 
non-relativistic matter and CMB with anisotropic pressure understood as 
the quadrupole moment of the distribution function for the photon gas fl^l 
- [^S) under a Kinetic Theory approach ]l9[ , p0| . For a matter-radiation 
mixture in which there is interaction between the components (a radiative 
photon-electron interaction), the anisotropic pressure can be interpreted 
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as shear viscosity within a causal irreversible thermodynamic approach || , 
pl| -|p2|. This provides an adequate model for the radiative era of cosmic 
evolution, from after nucleosynthesis to up to "matter" and "radiation" 
decoupling. 

In section II, Einstein's field equations are integrated up to a Friedmann- 
like equation, a quadrature in which the free parameters are three initial 
value functions related to the average of initial energy densities , pf 
and WR, the 3-dimensional Ricci scalar along the initial surface t = ti. In 
section III, we describe the possible physical interpretations of the models 
source focusing on the interacting matter-radiation mixture. In sections IV 
and V we express the free parameters in terms of suitable volume averages 
and initial contrast functions defined along t = ti, leading to simplified 
and elegant forms for all the relevant geometric and physical quantities. 
These quantities become fully determined (up to initial conditions) once 
the Friedmann-like quadrature is integrated (appendix A) yielding canoni- 
cal and parametric solutions. The models contain three important particu- 
lar cases presented in appendix B. In sections VI and VII, we examine the 
state variables of the model under small density contrasts. Theses small 
contrasts lead to two important classes of initial conditions equivalent to 
the definition of adiabatic and quasi-adiabatic perturbations of the initial 
densities and the initial curvature. In section VIII, we study the thermo- 
dynamic consistency of the coupled mixture model. 

Different variants of the imperfect fluid generalization of LTB solutions 
have been presented previously: using the equation of state of a non- 
relativistic ideal gas Jl4| , considering various possible ideal gas equations 
of state and a generalization to non-spherical geometry of the Szekeres- 
Szafron type |15 and the parabolic case for the matter-radiation mixture 
||l6|| . The present paper extend and complement the results of previous 
literature with respect to the specific cases of spherically symmetric and 
curved (elliptic and hyperbolic) models with a radiation-matter mixture 
source. 



II. FIELD EQUATIONS. 



Consider the usual LTB metric 

y/2 



1 -K 



dr 2 + Y 2 [d9 2 + sin 2 (6>)d</> 2 ] 



(1) 



where Y = Y(t, r), K = K(r) and a prime denotes derivative with respect 
to r. Instead of the usual dust source, we shall consider the stress-energy 
tensor of a fluid with anisotropic pressure 



T ab = pu a u b + ph ab + U ab 



h ab = c- 2 u a u b + , 



n a a = , 



(2) 



where the most general form for the anisotropic pressure tensor for the 
metric ([!]) with matter source (||) is given by: 11% = diag [0, — 2P, P, P], 
with P = P(t,r) to be determined by the field equations. Einstein's field 
equations for (Q) and (^) are 



8ttG 



Y(Y 2 + K) 



Y 2 Y' 



(3) 



~P = o( 2G 6 + G r) = 



Y(Y 2 + K) + 2Y 2 Y 



3Y 2 Y> 



(4) 
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Y(Y 2 + K) + 2Y 2 Y 

Fa 



(•5) 



where Y — u a Y a = Y :Ct . In order to integrate these equations we need 
to impose a relation between p and p but no previous assumption on P is 
necessary. 

Consider p in (|^) to be the sum of two contributions: non-relativistic 
matter described by dust plus radiation energy density (p^ ) 



p = P M 



P 



(r) 



1 (r) 
P = -/ \ 



with: 
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2M' 
Y 2 Y> 



(6) 



where M = M(r), while the form of p^ will be discussed further ahead. 
Inserting (@) into (ffl and ([!]) and integrating once with respect to t leads 
to the Friedmann equation 



Y 
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Y 



2M + W 



.Y 

Y 



K. 
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where W = W(r) and Yi — Y(ti,r) for an arbitrary fixed value t = U. The 
interpretation for M and W follows by substituting ([jj) into equations (^), 
(§ and (§) 
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w 
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Since Y 2 Y' is proportional to the determinant of the spatial part of the 
metric (Q), it is a covariant measure of proper local volumes. Therefore M' 
and W must have units of length and so it is convenient to define them in 
terms of the initial energy densities of matter and radiation 



2M = 



IttG 



» 



Y 2 Yjdr, W = 



IttG 



YfY'dr, 



(9) 



U 



where p^ n \ p\ are p^ m \ p^ evaluated at t 

For a non-rotating fluid with a geodesic 4- velocity the two nonzero kine- 
matic parameters are: the expansion scalar, O = u a . a , and the shear tensor, 
Cah = U( a ;b) ~ (Q/3)h a b- These parameters for (Q) take the form 

2Y Y' 



(10) 



a\ = diag[0,-2a,a,a], a = 1 ( ^ - ^ j , o ah o ah = 6a 2 , (11) 

Under the assumptions ^) the energy and momentum balance law, T ab ; & — 
0, associated with (|l|) and (|2|) gives 

p(r) + le P (r) + 6a P = 0, 

h b a {p, b + W b , c ) = => ( p -2P)'-6Py = 0, 

clearly illustrating how the pressure gradient is exactly balanced by the di- 
vergence of II a allowing for non-zero pressure gradients to be compatible 
with geodesic motion of comoving matter. 
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III. MATTER-RADIATION MODELS. 



The relation between p and p given by (^) suggests that the matter source 
can be understood as a mixture of non-relativistic matter (to be referred 
as "matter") and ultra-relativistic matter (to be referred as "radiation"), 
both characterized by the same 4-velocity. Therefore, a suitable physical 
interpretation for this source is an interactive, 'tightly coupled' mixture of 
matter and radiation examined under a hydrodynamical approach. The 
anisotropic pressure becomes a dissipative term (a shear viscosity) that 
must be examined within a thermodynamic framework pl| -p2|, p3| , ]2~j| . 
This type of source provides a convenient description for the "radiative 
era" , after nucleosynthesis and before decoupling of radiation and matter. 
An alternative approach is that of a decoupled mixture (non-relativistic 
matter plus the CMB) under the framework of Kinetic Theory jl^] - . In 
this paper we shall consider only the hydrodynamical approach applicable 
to the radiative era, leaving the study of a decoupled mixture for a future 
work. 

The 'tightly coupled' mixture of non-relativistic matter and radiation (a 
photon gas) is characterized by local thermal equilibrium (common tem- 
perature) among the components. This situation implies that interaction 
timescales are smaller than the cosmic expansion timescale, thus a hydro- 
dynamical approach is valid so that this interactive mixture behaves as a 
single dissipative fluid. Since heat flux necessarily vanishes for the LTB 
metric (0), this fluid can be described by the momentum-energy tensor 
(||) where U ab now becomes the shear viscous tensor. Considering that 
non-relativistic matter can be modeled by a classical monatomic ideal gas, 
while the photon gas satisfies the Stefan-Boltzmann law, the tightly cou- 
pled mixture of these components requires that p and p in (||) must comply 
with the equation of state 

3 1 
p = mc 2 n + -nkT + aT 4 , p = nkT + -aT 4 , (12) 

where k, a are Boltzmann and Stephan-Boltzmann constants, T is the com- 
mon mixture temperature, m is the mass of the most representative species 
of non-relativistic particles and n is particle number density, satisfying (we 
assume there is no net creation or annihilation of particles) the conservation 
law 

[nu a ), a = n = ^M., (13) 

where N(r) is an arbitrary function. If nkT -C aT 4 , or equivalently 

aT 3 /nk 3> 1, but mc 2 n/aT 4 is not negligible, then (12) can be approx- 
imated by [§,|5 

p = mc 2 n + aT A , p = \aT A , (14) 

Comparing the conservation law ([II]) and the equation of state (|l~4|) with 
(6a) and (6b), it is evident that the generalized LTB solutions can provide 
a model for the radiation-matter tightly coupled mixture if we identify 

p<«0 = mc 2 n, mc 2 N = 2 M' = ^± p^ Y 2 Y( , p\ m} = m<?n h 

(15) 

so that the radiation energy density and the temperature are given by 

= ^pir) = ^ aT 4 ^ aT ± = W'Y 2 Y' (16) 

Y 2 Y' c 4 c 4 c 4 
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As an alternative to (|12j) and (JL4|), it is possible to consider, instead of the 
Stefan-Boltzmann law, — aT 4 , the energy density and pressure of the 
radiation component as those of an ideal ultra-relativistic gas 

pM = nWfcr, P {r) = 3p (r) , « (r) = 775777 (17) 

where nA r ) is the corresponding particle number density, independently 
satisfying a conservation law like (la) with iV( r ) = N^ r \r). This approach 
has been followed previously in |^|and Jig] , the advantage being a simpler 
expression than (jig ) for the temperature in terms of the gradients of Y 



8irG , , 1 
— ^ Y = 3 



= N (r) kT, (18) 



In global thermal equilibrium, the fact that photon entropy per barion 
is conserved: — (4/3) aT 3 /n — const, implies that the ratio of pho- 
tons to baryons is constant and aT 4 cx nkT, hence the Stefan-Boltzmann 
and ideal gas laws are equivalent. Since we shall consider near equilib- 
rium conditions in which these two laws should be almost (but not exactly) 
equivalent, then aT 3 /(n k) is proportional to the (approximately constant) 
number of photons per non-relativistic particle. If the latter are baryons 
and electrons, then this quantity is very large, thus justifying the approxi- 
mation leading from ( p^ ) to (jlj). If temperatures are high enough so that 
the ratio aT 4 /(mc 2 n) is low enough for creation/annihilation processes 
to cancel each other (hence ( |l3"| ) holds), then ( |l4| ) provides a reasonable 
description of cosmic matter in the radiative era. The dominant radia- 
tive processes characteristic of this era mostly involve the photon-electron 
interaction (Thomson and Compton scattering, Brehmstrallung, etc) ||, 
PH and [p2| . The temperature range for the radiative era (up to matter- 
radiation decoupling) is roughly 4 x 10 3 A". < T < 10 6 A. The anisotropic 
pressure II a & is now the shear viscosity tensor, hence its evolution law must 
be consistent with the shear viscosity transport equation of a causal irre- 
versible thermodynamic theory. In the particular case when shear viscosity 
is the single dissipative flux, the entropy per particle and shear viscosity 
transport equation provided by Extended Irreversible Thermodynamics are 

s = s {e) + 4 F n ah n ah , => (snu a ). a = snu a > 0, (19) 
nl 



Tt\ cd h c a h d b + n ab 

where = (4/3)aT 3 / n is the equilibrium entropy per particle, taken (ap- 
proximately) as the initial photon entropy per non-relativistic particle, 77, r 
are the coefficient of shear viscosity and a relaxation time which, together 
with a, are phenomenological coefficients whose functional form depends 
on the fluid under consideration. A "truncated" version of (f20|), also known 
as the "Maxwell-Cattaneo" transport equation, is often used^2|, |23| for 
the sake of mathematical simplicity: 

Tt\ cd h c a h d b + H ab + 2 V a ab = 0, (21) 

Although this truncated equation satisfies causality and stability re- 
quirements, numerical and theoretical studies indicate that it might be 
problematic |2q|. In spite of the obvious limitation of having only one dis- 
sipative flux, the generalized LTB models provide arguments to infer the 
ranges of applicability of the truncated and full equations. 



2r]a ab = 0, 



(20) 
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The application of the thermodynamic formalism to the generalized LTB 
models, considered as a models of a hydrodynamical 'tightly coupled' mix- 
ture of matter and radiation, requires a convenient selection of the phe- 
nomenological quantities ry, a, r, and then solving the appropriate trans- 
port equation (either (pc| ) or (pl|)). In order to proceed with this task we 
will find in the following sections an intuitive characterization of initial 
conditions and suitable forms for the state variables. 



IV. INITIAL CONDITIONS AND VOLUME AVERAGES 



It is convenient to transform the Friedmann equation (Q) into the simpler 
quadrature 



-Ky 2 + 2fiy + lu 

y 



(22) 



where 



Y 



M 



y = v» ^ — v3 



YF 



W 

W 



K 

Y 

Y? ' 



Bearing in mind (j^), assuming that the initial hypersurface t = tj is ev- 
erywhere regular and contains at least a symmetry center 16 jL5|, we 
can express the quantities [A, u), k as the following type of volume averages 
along t = U 



(p 



fY*Y/dr 



2/j, 



8ttG (m) 



3c 4 



where 



(P 



/(3) 



(r)\ = fp^Y'dr 
~ fY?Y>dr 



Hi) = 



jW-RiYfYi'dr 



jYfY'dr 



2{KYi)' 
Y?Y< ' 



8ttG 



(P. 



^( (3) ^>,. 



(23) 



is the Ricci scalar of (1) evaluated at t = U and the range of integration goes 
from the symmetry center up to a comoving sphere marked by r. Notice 
that these averages (save for K = 0) do not coincide with the local proper 
volume defined by (|), namely: YfY'J ^/T~K = (l/3)d{Yf)/y / T^K, 
though they can also be characterized covariantly, since the metric func- 
tion Y (and so Y t ) is, in spherical symmetry, the "curvature radius", or the 
proper radius of the orbits of the rotation group SO(3). The volume aver- 
ages ( p3| ) lead to a compact and elegant mathematical description of initial 
conditions. From the volume averages in (pt) we can define the functions 

A. 



(m) A (r) A (fe) 



satisfying the following appealing relations 



P { l n) = (p^) 
pP = {ft* 



1 + A 



(m) 



1 + A 



(r) 



Wfti = (W^) [l + A 



(24) 
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justifying their interpretation as initial density and curvature "contrast 
functions" , as they provide a measure of the contrast of initial value func- 
tions , pf^ , ^TZi with respect to their volume averages along the initial 
hypersurface t = U. Because of their definition in terms of the volume aver- 
ages (|2^) , the value of the contrast functions at a given r — r* depend on the 
values of the functions p[ m \ p^, ^TZi on the integration range < r < r*, 
in which we assume that r = marks a symmetry center characterized by 
Y(t, 0) = Y(t, 0) = p| p|. Therefore, for a rest mass density lump p| m) 
decreases with increasing r, while for a rest mass density void it increases 
with increasing r, hence (from (pi|)) we have — 1 < A^ < for a lump 
and A| m) > for a void. The same criterion holds for the initial radiation 

density and scalar curvature, pf' and ^TZ-i. 

From ( ]25| ) and (|24|), it is straightforward to verify that the functions 
pt\p^\^Tli relate to A[ m) , A\ r \ A[ k) by 

w 2p[ r) 1 + A| m) k c 4 ^ 1 + A| m) 

£l .. IW1 . . It) ' £ 2 .. _ „ W . . Ik) ' V ^ 



A* pj m) 1 + Af M 87rGpl m) 1 + A| fe) 



and 



(r) w'/bJ ( m ) _ M7M .(fc) «'/« . . 

4 " YUYi 1 - Y!/Yi> ^ ~ Y'/Y- W 

Since the state variables P and kinematic parameters O, a are 

given in terms of quantities that will depend on the initial value functions 
/i, to, k and their gradients, it is useful to be able to express these gradients, 
by means of (^5|) and ( p6| ) , in terms of the more intuitively appealing initial 
contrast functions defined by (|24|). We examine this point in the following 
section. 



V. STATE VARIABLES AND KINEMATIC PARAMETERS 

Using (||) ( |l5| ) and ([!(]) , we can rewrite the state variables p^ , p^ , P 
in a more appealing form 

_( m ) 

(m) 2 Pi (m) 2 (c,n\ 

p y ' = mc n = , pi ' — mc m (27) 



y 3 r 



= 3p = ^-j^-, p\ r) = aT? (Stefan-Boltzmann), 



p^ = Zn^kTi (Ideal Gas), (28) 



T = — — (Stefan-Boltzmann), T = -— (Ideal Gas), (29) 

v V r / v 

p = w ' (30) 

where the auxiliary functions T, $, characterizing the spacial depen- 
dence of p, p, P, are given by 

r = Y 'I Y _ i , tfVy , T ,_ 4 + 3A^-r 4 + 3AW-4T 

W W 3[1 + Af»] ' 3[1 + A«] ' 

(31) 
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where has been defined in (]24|). The function T can be obtained as a 
function of y, fi, u>, n and the initial contrast functions by performing the 
integral quadrature of (f22[). Let 



z ^ , ydy = — [ ydy (32) 

' ' - K y 2 + 2ny + Lo} 1/2 [-e 2 y 2 + 2 y + £l ] 1/2 



where 



Zi(H,U,f) = Z \y=l, e l = -j £ 2 = -, 

M A 4 

then, by using the chain rule on Z — Zi we can obtain for T, 4" and $ in 
( |3l| ) the following expressions 

r = 1 - 3AAf° - 3BAf - ZCAf\ (33) 



1 + AA| ro ' + (l + flJA^ + CA^ 
* = ! j[ + A (D ( 34 ) 



$ = £AAg + (l+jg)A^ + 4CAf 
1 + A 2 W 

where 

. n[d(z - Zi)/dn] „ uj{d(z~ Zi)/du>] K [d(z- Zi)/dK] 

A = r „ — ,„ : , is = r „ — , - : , O = 



y[dZ/dy] ' " y[dZ/dy] ' 

(36) 

In order to obtain the explicit functional forms for A,B,C above we need 
to evaluate Z explicitly This is done in Appendix A, while a convenient 
interpretation as "initial contrast functions" has been provided in the pre- 
vious section for the functions A^ , A^ , A^ . It is important to mention 
that no a priori assumption on P was made in order to obtain (|30|). This 
specific form of anisotropic pressure (shear viscosity) follows directly from 
g, (§, © and (§). 

In terms of the new variables, the kinematic parameters, 0, a, introduced 
in @ and (|ll|), take the form 

3y t 
© = — + p 

y r 

= VQ 1-(3A + yA.^A^ - (3B + yB, y )Aj r ] - (3C + yC, y )A^ 
y 1 - 3AA„ (m) - 3BA, (r) - 3CA„ (fc) 



CT £_ _ VQ_ A^+B^+CyA^ 

3r y 1 - 3AA, (m) - 3BAf } - 3CAf } ' 

where Q = — ny 2 + 2jiy + LO and a sub-index y means partial derivative with 
respect to y. The new solutions presented so far become fully determined 
once the Friedmann equation ( p2| ) is integrated (Z is explicitly known) 
for specific initial conditions provided by fi, lu and k. This integration is 
presented in appendix A and particular cases of interest are presented in 
appendix B. The conditions for the models to comply with regularity at 
the center and energy conditions are given in |[6[ and p"5| . 
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VI. SMALL DENSITY CONTRASTS AND AN ADIABATIC 

EVOLUTION. 

The thermodynamic study of the generalized LTB models, under the 
framework of Extended Irreversible Thermodynamics, necessarily assumes 
small deviations from equilibrium. Since the equilibrium limit of the models 
(excluding the dust limit which has trivial thermodynamics) is the FLRW 
limit that follows by setting the initial contrast functions to zero, then a 
"near equilibrium" evolution can be related to a "near FLRW" evolution 
that follows by assuming "near homogeneous initial conditions" , or in other 
words, small initial contrast functions, ie: 



implying 



a^i « i, i An « i, i An « i, 



(39) 



2/.t 



SttG ( m ) 
3c 4 Pl ' 



8itG r r ) 



.(3) 



Hi, 



2p 



(r) 



Am) 



(2 



MTU 



8irG J m ) 

ri 



(ei,£2 



const. 



(40) 



where the exact form of £1,62 is given by (|25j). Under these assumptions, 
all functions depending on (n,u),K,y) become approximately functions of 
y only. From §, @, M), © and pl), this implies 



Z « Z(y), t-Un Z{y) - Zi, A « A(y), B » B(y), C « C(y), (41) 

with Zi rs const., so that hypersurfaces y — const, approximate hypersur- 
faces t — const, and y becomes the time parameter. The radial depen- 
dence in quantities like T, 'J, <£> is then contained in the initial contrast 
functions. Regarding the state variables, a Taylor series expansion around 
A M A M A (fe) yields at first order 



r 



4 . $ 

3 5 ' r 



A 



(r) 



where we can formally identify the "perturbation" as 



5 = 3[AA( m) +BA| r) 
so that, from ( p7| ) and (|4^), we have 



CAf»] = i-r, 



(42) 



(43) 



[1 + S], P 



(r) ~ 



M r 



P 



(r) 

6y 4 



A 



(r) 



(44 



Therefore, under the "small contrast approximation" given by (|39|), ( j40| ) 
and ( fidf ) the models follow a nearly homogeneous evolution with respect 
to a scale factor y « y(t), deviating from homogeneity in a way (as given 
by equations (0)) that is formally analogous to that of adiabatic pertur- 
bations. This situation holds for whatever form of the contrast functions, 
as long as these adimensional quantities are very small (so that the initial 

functions P { [\ are almost constant). Notice that, at first order 
expansion on the contrast functions, we have for (E9j) 



1/4 



T 

y 



(45) 
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so that both, the Stefan-Boltzmann and Ideal Gas laws, yield the same 
expression for T. Therefore, as long as we assume the small contrasts 
approximation, the Stefan-Boltzmann and Ideal Gas laws yield the same 
result and we can use them indistinctly. 

The resemblance of to expressions characteristic of adiabatic pertur- 
bations is not surprising: the generalized LTB models under consideration 
(under a thermodynamic approach) have zero heat flux and so can be as- 
sociated with thermodynamic processes that are indeed adiabatic, though 
irreversible (because of the shear viscosity). The connection between the 
contrast functions and adiabatic perturbations on a radiation-matter 'tight 
coupling' is an interesting feature that deserves proper examination and 
will be studied in a separate paper. We will use this analogy as a theoret- 
ical tool for understanding the type of quasi-homogeneous evolution and 
for defining suitable initial conditions for the study of the thermodynamic 
consistency of the models. 



VII. ADIABATIC AND QUASI-ADIABATIC INITIAL 
CONDITIONS. 

As shown in p6j, t he quantity Aj*' = A^ — (4/3) A^™ 1 ' becomes, under 
the assumptions (39), ([h]) and (|4^), a sort of average change of photon 
entropy per barion at the initial hypersurface. Hence, it is convenient to 
rephrase initial conditions in terms of this quantity. After some algebraic 
manipulation on the explicit general forms given in Appendix A for the 
functions A, B, C defined by (^6|), we find that the following restrictions on 
the initial contrast functions 



3 1 



At) 



A 



A 



(fe) 



1 



» 



2^ ■ 



(46) 



leads to the following compact forms for the function T defined in (33) 

V? 1 A (r) 



r = i 



l - 



iej + Ae lV - y 2 



3 

4y2 

for the parabolic case (ea = 0), 



A| (8e? + 4ei -l) 
v it 



,(») 



y 2 Vql 



A 



(r) 



4e 2 A 2/ 2 



ei 1 



'^'■'1.-41 



A (r? - rii)y/q 



A 



(<0 



elliptic (e 2 > 0, +sign) and hyperbolic (e 2 < 0, — sign) cases, (47) 
where e\ = u>//j,, e 2 = «//-t where defined in ( p5| ) and 

q = -e 2 y 2 + 2y + ei, qi =-e 2 + 2 + ei, 
Aq = l±ei|e 2 |, elliptic (+ sign), hyperbolic (— sign), 



7] = arccos 



1 - £2 2/ 



77 = arccosh 



, r/i — arccos 



l-e 2 



rji = arccosh 



i + N 



, elliptic case, 
hyperbolic case, 
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It is evident, by looking at these forms for T above, that this function takes 
a very simple form 



1 



1 - 



A 



(r) 



(48) 



valid for all cases (parabolic, elliptic and hyperbolic) if we use initial con- 
ditions given by 



A 







A. 



(m) 



3 A (r) 

4 1 ' 



A [k) = -A 
2 



1 



(r) 



(49) 



Since all state variables, kinematic parameters and auxiliary functions arc 
constructed from T, y and initial value functions, the assumption ( [49| ) 
leads to very simplified forms for all expressions. Considering the initial 
contrast functions A[ m \ A,[ r \ A^' in ( p^ ) as formally analogous to exact 
perturbations on p\ m \ p| r \ < - 3 - ) Tli, the factor 4/3 relating A^ and Aj in 
( [49| ) is reminiscent of adiabatic perturbations on the initial value functions. 
Hence, following Jll| and pq] , we will denote ([l9]) as "adiabatic initial 
conditions" , while the more general case (f|| with A[ s) ^ will be refered 
to as "quasi-adiabatic initial conditions" . For the remaining of this paper 
we will only consider initial conditions of either these two types, under the 
"small contrast approximation" given by (^), (^0|) and ([4l|). 



VIII. THERMODYNAMIC CONSISTENCY 



In order for the models presented in this paper to be physically meaning- 
ful they must satisfy energy conditions and must be compatible with causal 
Extended Irreversible Thermodynamics. Before to start this analysis we 
need first to provide an expression for the phcnomenological quantities 
77, a, t (coefficient of shear viscosity, relaxation time). Following the ap- 
proach used in |l6| and p6|] , we consider for a matter-radiation mixture 
interacting via radiative processes, the "radiative gas" model, associated 
with the photon-electron interaction, which provides the following forms 
for 7} and a |H| -|| 



4 4 4 

„ = - pT = -aTr, 



T 



15 



8aT 4 ' 



Inserting (BOT) into (|19|) and rt2Q), the latter equations become 



5 k U ab U ab 
8r p 2 



15 k P 2 
At p 2 '' 



(50) 



(51) 



P 



P 

vq 1 — 

,P 



0, 



(52) 



where in ( |52| ) v$ = 5/4, 25/32, respectively, for the ideal gas and Stefan- 
Boltzmann laws (we examine the truncated equation ( |2l| ) further ahead). 
The transport equation (^2|) is an evolution equation for P imposed by 
thermodynamic theories external (even if coupled) to General Relativity. 
On the other hand, there are evolution equations for P that follow from 
the field equations, for example, by evaluating P from (||) and eliminating 
Y/Y and Y'/Y' with the help of @ and ([ll]), leading to 



P 



5a+-Q) P - 2pa = 0, 



(53) 
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an exact equation that must be satisfied by the generalized LTB models. 
Obviously, an equation like ( |5^ ) does not exactly coincide with (|5^), and 
so compatibility between these models and Extended Irreversible Ther- 
modynamics requires finding an appropriate expression for r that should 
make (|5^) consistent with the transport equation. Comparing ( |52| ) (with 
Vq = 5/4, ideal gas) with the evolution equation (|53"1), we can see that both 
equations coincide if we identify 

4 + 3A - r) - r 



a {2P 2 + fp 2 -5pP} 4er 4 

This expression is justified as long as it behaves as a relaxation parameter 
for the interactive matter-radiation mixture in the theoretical framework 
of EIT as we will explain bellow. 

It is also useful to compute the collision times for Thomson and Compton 
scattering (the dominant radiative processes in the radiative era) 



-pP 



3A 



4r 



4 + 3A« + §r 



~ 256 



t = 



m P c 
t c = -r^—t 



(55) 



where ut is the Thomson scattering cross section, m e is the electron mass 
and n e is the number density of free electrons, a quantity obtained from 
Saha's equation, leading to 



1 

7 ~ 2ca rr n( m ) 



4/iV m )exp (B /k g T) 



1/2 



(2-7T m e k B T) 



3/2 



(56) 



where Bq and h are respectively the hydrogen atom binding energy and 
Planck's constant. For details of the derivation of (54) see flj6| . 

The restrictions for physical acceptability and thermodynamic consis- 
tency were discussed in [^6| and fill . We provide a summary in the follow- 
ing list: 

• (i) Positive definiteness and monotonicity of and p^ m \ and \P\ <C 
V- 

• (ii) Regularity condition 

T > 0, (57) 
which prevents the occurrence of a shell crossing singularity [jljj . 

• (iii) Positive definiteness of r and s (consistent with positive entropy 
production). The sign of s (see ([si])) depends only on the sign of r. 

• (iv) Concavity and stability of s by requiring 



t > 0, T 

s 



2<jT \P- 



3*$ n ( r ) 

ri 



PI 



3p 



(r) 



^<0. 
T 



(58) 



• (v) Appropriate behavior of the relaxation time r. 

During the radiative era the Thompson time £ 7 must be smaller than 
the expansion (or Hubble) time, approximately defined by tu — 3/0. 
The decoupling hypersurface is then defined by i 7 = tu ■ From (iii) 
and (iv) above jL5|, Jl6|, t must be a positive and monotonously 
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increasing function (if the fluid expands), then it must also be (dur- 
ing the interactive period) qualitatively similar but larger than the 
microscopic timescales of the various photon-electron interactions oc- 
curring in the radiative era. It is usually assumed that r is of the 
order of magnitude of a collision time, and as matter and radiation 
decouple all these timescales must overtake the Hubble expansion 
time t H = 3/0. A physically reasonable r should be comparable in 
magnitude to t 7 near decoupling and must have an analogous qual- 
itative behavior to t 1 . Therefore r must be smaller than tn during 
the radiative era and then must overtake it at the decoupling hyper- 
surface. 

In the following list we test the previous restrictions (i) to (v) under 
the assumption of the small contrasts approximation (|39|) , ( p0| ) and , 
together with either adiabatic ( f49| ) or quasi-adiabatic ([46|) initial conditions: 

• (i) For a wide range of initial conditions all these quantities have a 
physically meaningful behavior as indicated in Figure. [I]. 

• (ii) Figure. |^ shows that V is positive and almost equal to unity for 
most of the evolution range of y. For the small range defined by 
1(T 5 < |A^ S) | < 1(T 3 and log 1Q (j/) < 1 the plot of T is distinct for 
lumps and voids. 

• (hi) and (iv). Figure. || shows that r > 0, s > are satisfied for 
quasi-adiabatic initial conditions. From Figure. [3] it can also be seen 
that r is a monotonically increasing function and so s < 0. The same 
results hold for adiabatic conditions (see Figure, ^a). 

• (v) Under adiabatic initial conditions (|A^| = 0), r does not have 
the required behavior described above (never overtakes t H ). This is 
shown in Figure. ||a. We find that the desired behavior of r is encoun- 
tered for quasi-adiabatic initial conditions. As shown in Figure. ||, this 
happens in the elliptic case for the values 1CT 55 < |A^| < 10^ 3 \ 
compared to |A^ S ^| < 10^ 8 ^ in the hyperbolic (not shown) and 
parabolic case ^6|. Finally, we plot the ratios r/t H , t c /t H and t 1 /t H 
in Figure. ^| showing that these times have a physically reasonable 
and consistent behavior. 

Furthermore, we examine the truncated transport equation for the gen- 
eralized LTB models that follows by inserting (J50|) into (pT 



P+ -P+ -pa = 0, (59) 
t 5 



A comparison between fl52| ) and (|59j) shows that the full equation is equiva- 
lent to the truncated one if P/p -C 1 and 1/r <€. 49/3. The first condition 
simply requires small deviations from equilibrium (as shown by (|5l|)) and is 
compatible with energy conditions [^5). However, the condition 4r 3> 3/0 
is more problematic and is only reasonable after matter and radiation have 
decoupled, hence, under the assumptions underlying the thermodynamic 
study of the models, the truncated equation cannot describe the interact- 
ing period nor the decoupling process. This shortcoming of (Elf) has been 



discussed previously for the parabolic case 26 . As shown by Figure. |6p, 
this situation holds also for the elliptic and hyperbolic cases. 

In order to complete the analysis we compute the Jeans mass associated 
to the initial conditions of the case A}' 7^ 0. This mass is given by|9j,||,[0] 
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(60) 

, Xl = PT/Ph * and r 
follow from (|33|) and C s is the speed of sound, which for the equation of 
state (|l4|), has the form 

-l 

, p< m > = mcV m \ p^ = 3n^k B T. (61) 

Evaluating © for y = y D w 10 2 4 , e w l/x< « 10 3 and p^ w a B T/ w 
7.5 x 10 9 ergs/cm 3 , yields Mj ss 10 49 gm, or approximately 10 16 M©. This 
value coincides with the Jeans mass obtained for perturbative models domi- 
nated by baryons in the radiative era as decoupling is approached. Finally, 
it is worthwhile remarking that we are dealing with an epoch where a 
cosmological constant or a quintessence component should be clearly sub- 
dominant |27j and ignored. 

IX. CONCLUSION 

We have derived and discussed important generic properties of a class 
of exact solutions of Einstein's equations that generalize the famous LTB 
solutions with a dust source to an imperfect fluid with anisotropic pressure 
H a b- The integration of the field equations does not involve Tl a b, though 
once this integration is done this pressure becomes also determined (up to 
initial conditions). The issue regarding the compatibility of H a b and its 
evolution law with the physical assumptions underlining the models have 
been addressed. These models provide a physically plausible hydrodynam- 
ical description of cosmological matter-radiation mixture in the radiative 
era, between nucleosynthesis and decoupling. By assuming small initial 
density contrasts (consistent with small deviations from equilibrium), we 
have shown that the state variables of the models are qualitatively and 
formally analogous to that of adiabatic perturbations on a FLRW back- 
ground. We have also found two classes of initial conditions, based on well 
defined initial contrast functions, that are formally equivalent to the defi- 
nition of adiabatic and quasi-adiabatic perturbations on initial value func- 
tions p\ m \ p\ r \ In particular, we showed that for quasi-adiabatic 
initial conditions the models are thermodynamically consistent and physi- 
cally acceptable. However, this consistency does not hold for the truncated 
transport equation (as shown also in pq|), only for the full transport equa- 
tion of Extended Irreversible Thermodynamics the relaxation time of shear 
viscosity has the appropriate physical behavior. 
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APPENDIX A: INTEGRATION OF THE FRIEDMANN 
EQUATION (22). 

The quadrature (^ ) depends on the sign of n = KjYf. The integral 
yields either an implicit function ct = Z{y, fi, oj, n) where Z is defined by 
( |32|) (the "canonical solution"), or a parametric solution having the form 
y = yjji, r), t = t{rj,r). The canonical solution of the equation equivalent 
to (|2|) found in the literature for the LTB dust case |l| , g] , is usually given 
in the form c[t — to(r)] = Z(y,fi, k), where the function to(r) comes as an 
integration constant. Instead, we will provide the canonical solutions in 
the form c{t — tj) = Z — Zi, so that we can identify to(r) = U + Zi. The 
function to (r) can also be identified as the "big bang time" . Only for k = 
the canonical solution can be inverted as y — y(t, r) by solving a cubic 
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equation, however the resulting expression is too cumbersome and will not 
be given. 

All variables introduced in previous sections will be made fully deter- 
mined as functions of y and of r, with the dependence on r given through 
the initial value functions and their initial contrast functions. As usual in 
the study of dust LTB solutions, we will classify the integrals (|3^) in three 
cases: "parabolic" (k — 0), "elliptic" (k > 0) and "hyperbolic" (n < 0). 
We shall examine each case separately below. 



1. Parabolic case. 

If k — 0, so that ^TZi — = 0, then ( p2[ ) leads to the canonical form 
4 



c(t -ti) = Z - Zi 



s]2y + e t (y - Cl ) - (1 - e x ) (Al) 



where e\ was defined in (p5|). The solution in terms of the parameter rj is 
1 



V = 2 



V 

ei 



c{t-U) = ^ 
bp 



V 

ei 

A* 



The functions A, B in ( p6[ ) have the form 

i ry^+r" 



.4 



3y 2 



V2T7 



i(l-4 ei -8e 2 )-(y 2 -4 ei y-8e 2 ) 



(A2) 



B 



Notice that the limit p^f' — > (or ei — > 0) leads to the known parabolic 
solutions of the dust LTB solutions, though, the limit p^ — > is singular. 



2ei 

y 2 



n . ^ \/2y + e! 

( 1 + e v—^=^=- - (y + e i) 



(A3) 



2. Elliptic case 
For k > the canonical integral (|32|) is 



c(i - = Z - Zi 



(A5) 



Q = — «y 2 + 2py + ui, Q :V = —2ny + 2/i, A = i(p 2 + kui) 
while the parametric solution is given by 



V : 



1 ^ 

1-— cos, 



(A6a) 



c(t - ti) 



,3/2 



The functions A, B, Cm (p&j become 



(A6b) 



A 



K U 



(Vi-V)VQ 4p 2 + A + Apu y/Q {Ap 2 + X)y + Apuj 



A 



/Q7 



A 



(A7a) 
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B = 



2uj 



[Alb) 



r ,_ SfiVQi^-U) [(A-4 / ^) K -(2 M + o ) )(4/i 2 + 2A)] 

^ — ~ — — ~ 1 



2k 3 / 2 ?/ 2 



2kAj/ 2 VQI 



[(A - ^Lu)ny - (2fj,y + w)(4/z 2 + 2A)] y/Ql 



(A7c) 



rj = arccos 



As with the parabolic case, we obtain the dust LTB elliptic solution in the 
limit lo —> 0. It is also interesting to compare the evolution of y as function 
of X] with the dust case. This evolution is also time-symmetric, but it begins 
at the value t]bb = arccos (2/i/VA) > 0, instead of t\bb — 0. The maximal 
value of y as the fluid bounces is y ma x = (/V K )(1 + larger than 

the value for the dust case y ma x — 



3. Hyperbolic case 

For k < the canonical integral ([32]) is 



c(t - ti) = Z - Zi = 



UI3/2 



arccosh 



, A > (A8a) 



c(t-ti) = Z-Z, 



VQ m ln f Q,» 

N M 3 / 2 Uv^k 



, A < (A8b) 



Q = |k|j/ 2 + 2it?y + w, Q >a = 2|k|j/ + 2/U, A = 4^ 2 - 4\k\oj 
while the parametric solution is given by 



^ v. i 
cosh 71 — 1 

2/1 ' 



(A9a) 



c(t - fi) 



13/2 



v^A . 

sinn ri — n 

2/i 



(A9b) 



The functions A, £?, C in (14) become 



A = 



(Vi-V)VQ , 4(i 2 + A + 4^w y/Q (4^ 2 + X)y + i^uj 



B 



2lu 

V 



(w + /i) 



Q 



/ft 



A 



(AlOa) 
(A106) 



^_ ^VQim-v) . [(A-4^) K -(2 M + W )(4 M 2 + 2A)] 
c^ — ~ — rrTo — - 1 



2|k| 3 /V 



2kAj/ 2 
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[(A - Afuojny - (2py + uj)(4p 2 + 2A)] 



(Aide) 



/; = arccosh j — Q= 



As with the parabolic case, we obtain the dust LTB elliptic solution in the 
limit cij — > 0. 



APPENDIX B: PARTICULAR CASES. 



Expressing initial value functions in terms of the initial averages and 
contrast functions defined in sections III and IV greatly facilitates the de- 
scription of the correspondence to the following particular cases: 

• Dust limit. If = {p { p) = 0, then w = and Q = (-Ky + 2p)y, 

hence we have: p = p^ = p["^ /(y 3 T) and p = P = 0, and so the 
solutions reduce to the "usual" LTB solutions with a dust source. 
The initial density and curvature contrasts A| m ' ) , have the same 
interpretation as in the general case. Notice that B — (from (14a)) 
and that the auxiliary functions 'J, $ are irrelevant in this sub-case. 

• FLRW limit. Homogeneity along the initial hypersurface follows 
by demanding A l (m) = A^ = Af } = 0, leading to p| m) = (^ ro) ), 

p[ r) = (pfV = ( {3) TZ,), so that p\ m \ p ( [\ ^Tli are constants. 
Equation (E2h implies then that p, uj, k are also constants, so that 
M, W (x Yfand K cx Yf. Therefore y obtained by integrating ( p2| ) 
must be a function of t only and Y — R(t)Yi. From ( |3^ ) to ( |35| ) we 
have: V = * = 1 and $ = leading to p = p(t), p = p(t), P = 
and a = 0, = 3y/y = 3R/R. Thus the particular case in which 
initial density and curvature contrasts vanish is the FLRW limit of 
the solutions, a FLRW cosmology where p, p satisfy the "dust plus 
radiation" relation (^|). 

• Vaidya limit. If p^ 1 ' > 0, p+ > for < r < r& but vanish for 

r > Tb, then (/O-™^) and (p^), as well as M and W are constants 
for r > ri, and we have (along this range) a sub-case of the Vaidya 
metric characterized by the mass function m — M(r&) + W{rb)Yi/Y . 
This particular Vaidya solution is the "exterior" field for the models 
derived in this paper, generalizing the Schwarzschild exterior of LTB 
dust solutions. Notice that as the fluid expands the Vaidya mass 
function, m, tends to a constant Schwarzschild mass given by Af(r&). 
More details are provided in Jl5| . 

The solutions derived in this paper can be smoothly matched along a 
comoving hypersurface marked by r — rb with their FLRW and Vaydia sub- 
cases. The matching conditions are discussed in detail in pl|. Just as LTB 
dust solutions can be generalized to the Szekeres dust solutions without 
isometries, LTB solutions with an imperfect fluid source can be gen eralized 
to the Szafron-Szekeres metrics with an imperfect fluid source |l5| . Even 
without spherical symmetry, these solutions can still be smoothly matched 
to the FLRW and Vaydia sub-cases along a spherical comoving boundary 
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FIG. 1: (a) Logarithmic plot of the densities ratio. These densities are positive 
definite and decreasing functions of log 10 (j/); The plot also shows that for the early 
stage the radiation density dominates while for later stages the matter density 
dominates, (b) The ratio P/p is very small complying with energy condition. P 
is negative/positive for lumps/voids respectively. As we are dealing with small 
initial density contrasts defined by ([w]), ( fio| ) and (El]), we used for all the plots 
ex « 2 P t ) / P \ m) = 10 3 and e 2 » (c 4 /8^G)( (3) %/p< m ') = 10" 7 . 



FIG. 2: The function Y given by ( p3[ ) in terms of log 10 (j/) and log 10 (|A^|), from 
the initial hypersurface y = 1 corresponding to Tj = 10 6 Kelvin. F is positive as 
required, it is almost equals to one for all the evolution range, (a) The lumps 
case, (b) the voids case. 



FIG. 4: Plot of log 10 (t/*h) in terms of log 10 (y) and log 10 (|A< s) |); (a) for 1CT 8 < 
|A' S ' < 1CF 5 ' 5 ) the relaxation time is not initially smaller than the Hubble time ; 
(b) For the range 10 -5 ' 5 < |A' S '| < 10'~ 3 ' the relaxation time is initially smaller 
than the Hubble time, and it overtakes it at about y = 10 2 ' 4 corresponding to 
Td = 4000/i' satisfying the required behavior. 



FIG. 5: The figure depicts the ratios iog 10 (t c /t H ), iog w (r/t H ) and log 10 (t-,/t H ) 
versus log 10 (y) for A' r ' = 10~ 4 and A' s ' = 10 -3 . Initially, the ratio r/t H is far 
larger than t 1 jt H but become comparable to it near the decoupling hypersurface 
defined by t 7 = t H . Note the similarity in the qualitative behavior of these 
two ratios. This plot shows also how the Compton time is important for high 
temperatures (near log 10 (y) = 0) but rapidly overtakes t H and so Compton 
scattering is no longer an efficient radiative process. Thompson scattering is very 
small initially but becomes the dominant process near decoupling. The relaxation 
time r is a mesoscopic quantity that acts roughly as an "average" timescale for 
all these processes. 



FIG. 6: Plot of \og 10 (r/t H ) in terms of log 10 (j/), (a) for adiabatic conditions 
|A' S '| = 0, r/t H is always smaller then 1, so r does not overtake t H , (b) For r 
defined from the truncated equation (49c) r/t H is always smaller then 1 as well. 



FIG. 3: Plot of log 10 (r) in terms of log 10 (j/) and log 10 (|A' |), r is a positive 
definite and increasing function for 10 -8 < |A' 3 '| < 10~ 3 and for all evolution 
range. 
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